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Abstract 

Modern techniques of the renormalization group (RG) combined with effective field theory (EFT) methods are 
revolutionizing nuclear many-body physics. In these lectures we will explore the motivation for RG in low-energy 
nuclear systems and its implementation in systems ranging from the deuteron to neutron stars, both formally and in 
practice. Flow equation approaches applied to Hamiltonians both in free space and in the medium will be empha- 
sized. This is a conceptually simple technique to transform interactions to more perturbative and universal forms. An 
unavoidable complication for nuclear systems from both the EFT and flow equation perspective is the need to treat 
many-body forces and operators, so we will consider these aspects in some detail. We'll finish with a survey of current 
developments and open problems in nuclear RG. 

Keywords: Renormalization group, nuclear structure, three-body forces 



1. Overview 

The topic of these lectures is the use of renormal- 
ization group (RG) methods in low-energy nuclear sys- 
tems, which include the full range of atomic nuclei as 
well as astrophysical systems such as neutron stars. We 
will examine why the RG has become an increasingly 
useful tool for nuclear physics theory over the last ten 
years and consider how to apply RG technology both 
formally and in practice. Of particular emphasis will be 
flow equation approaches applied to Hamiltonians both 
in free space and in the medium, which are an accessible 
but powerful method to make nuclear physics computa- 
tionally more like quantum chemistry. We will see how 
interactions are evolved to increasingly universal form 
and become more amenable to perturbative methods. A 
key element in nuclear systems is the role of many-body 
forces and operators; dealing with their evolution is an 
important on-going challenge. 

The expected background for these lectures is a thor- 
ough knowledge of nonrelativistic quantum mechanics, 
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including scattering, the basics of quantum field the- 
ory, and linear algebra (it's all matrices!). We will not 
assume a knowledge of nuclear structure or reactions, 
or even many-body physics beyond Hartree-Fock. No 
advanced computing experience is assumed (although 
Mathematica or MATLAB knowledge will be very help- 
ful in exploring simple RG examples). 

By necessity, we will only scratch the surface in these 
lectures. For a thorough treatment of flow equations for 
many-body systems not including nuclei, see the book 
by Kehrein (Q. For more details on applications of 
flow-equation and similar renormalization group meth- 
ods to low-energy nuclear physics, the review article 
and references therein are recommended. 



2. Atomic nuclei at low resolution via RG 

2. 1. Goals and scope of low-energy nuclear physics 

The playing field for low-energy nuclear physics is 
the table of the nuclides, shown in Fig.fT] There are sev- 
eral hundred stable nuclei (black squares) but also sev- 
eral thousand unstable nuclei are known through exper- 
imental measurements. The total number of nuclides is 
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still unknown (see the region marked "terra incognita"), 
with theoretical estimates suggesting it could be as high 
as ten thousand! These unstable nuclei are the object 
of scrutiny for new and planned experimental facilities 
around the world. The challenge for low-energy nuclear 
theory is to describe their structure and reactions. We'll 
return in the final lecture to discuss the overlapping re- 
gions where the theoretical many-body methods listed 
in the figure can be applied. 

Let's start with some of the questions that drive low- 
energy nuclear physics research. These include general 
questions about the physics of nuclei |3|: 

• How do protons and neutrons make stable nuclei 
and rare isotopes? Where are the limits? 

• What is the equation of state of nucleonic matter? 

• What is the origin of simple patterns observed in 
complex nuclei? 

• How do we describe fission, fusion, and other nu- 
clear reactions? 

These topics inform and are in turn illuminated by ap- 
plications to other fields, such as astrophysics, where 
one can ask: 

• How did the elements from iron to uranium origi- 
nate? 

• How do stars explode? 

• What is the nature of neutron star matter? 
or of fundamental symmetries: 

• Why is there now more matter than antimatter in 
the universe? 

• What is the nature of the neutrinos, what are their 
masses, and how have they shaped the evolution of 
the universe? 

Finally, there are applications, for which we are led to 
ask: How can our knowledge of nuclei and our abil- 
ity to produce them benefit humankind? The impact is 
very broad, encompassing the Life Sciences, Material 
Sciences, Nuclear Energy, and National Security. 

In Figure [2j the energy scales of nuclear physics are 
illustrated. There is an extended hierarchy, which is a 
challenge, but also an opportunity to make use of ef- 
fective field theory (EFT) and renormalization group 
(RG) techniques. The ratio of scales can become an ex- 
pansion parameter, leading to a systematic treatment at 
lower energies. The progression from top to bottom can 




Figure 1: The nuclear landscape. A nuclide is specified by the number 
of protons and neutrons |3 1. 
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Figure 2: Hierarchy of nuclear degrees of freedom and associated en- 
ergy scales |3|. 
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Figure 4: Momentum space representation of the Argonne vis (AV18) 
NN potential in the ^o channel 0. 



Figure 3: Some phenomenological potentials that accurately describe 
proton-neutron scattering up to laboratory energies of 300 MeV |4|. 



be viewed as a reduction in resolution. In these lectures, 
our focus is on the intermediate region only, where pro- 
tons and neutrons are the relevant degrees of freedom. 
But even within this limited scope, the concept of reduc- 
ing resolution by RG methods is extremely powerful. 

2.2. Lowering the resolution with RG 

What do we mean by resolution? Even the general 
public these days is familiar with the concept of dig- 
ital resolution for computer screens, cellphones, tele- 
visions. High resolution is associated with more pix- 
els, which follows because as pixel size becomes small 
compared to characteristic scales in an image, greater 
detail is seen. In our discussion, we associate resolution 
with the Fourier transform space and the phenomenon 
of diffraction. 

Recall the basic physics: if the wavelength of light is 
comparable to or larger than an aperture, then diffraction 
is significant. If there are two sources, we say we can 
resolve them if the diffracted images don't overlap too 
much. For a fixed angle between sources or details in 
the object being observed, we find that the wavelength 
determines whether or not we resolve the details. Being 
unable to resolve details at long wavelength is generally 
considered to be a disadvantage (e.g., for astronomical 
observations), but we turn it to an advantage. 

A fundamental principle of any effective low-energy 
description (not restricted to nuclear physics!) is that 
if a system is probed at low energies, fine details are 



not resolved, and one can instead use low-energy vari- 
ables for low-energy processes. Renormalization theory 
tells us that the short-distance structure can be replaced 
by something simpler without distorting low-energy ob- 
servables. The familiar analog from classical electro- 
dynamics is the replacement of a complicated charge 
or current distribution with a truncated multipole ex- 
pansion. In the quantum case, the replacement can be 
done by constructing a model, or in a systematic way 
using effective field theory. We emphasize that while 
observable quantities (such as cross sections) do not 
change, the physics interpretation can (and generally 
does) change with resolution. What if there is no exter- 
nal probe? Then the particles still probe each other with 
resolution set by their de Broglie wavelengths. Low- 
density nuclear systems would seem to imply low reso- 
lution. But the picture is complicated by the nature of 
traditional internucleon potentials. 

Figure [3] shows several phenomenological potentials 
that reproduce nucleon-nucleon scattering phase shifts 
up to about 300 MeV in lab energy. They are charac- 
terized by a long-range attractive tail from one-pion ex- 
change, intermediate attraction, and a strongly repulsive 
short-range "core". For our purposes, the partial- wave 
momentum-space representation, such as 

(k\V L =oW) ex Jr 2 drj (kr) V(r) j (k'r) (1) 

for S-waves is more useful. Here k and k f are the rel- 
ative momenta of the two nucleons. This is shown for 
the AVI 8 potential in Fig. in the l S channel. (The 
spin and isospin dependence of the nuclear interactions 
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Figure 5: Alternative momentum space representation of the AVI 8 
potential in the ! So channel |2|. 



is very important, unlike the situation in quantum chem- 
istry.) In the momentum basis, the potentials in Fig. [3] 
are no longer diagonal, so we need three-dimensional 
information, but we generally use the flat contour repre- 
sentation of the same information, as in Fig. [5] The RG 
evolution of potentials will be visualized as changes in 
such pictures. 

We work in units for which h = c = 1. Then the 
typical relative momentum in the Fermi sea of any large 
nucleus is of order 1 fm -1 or 200 Me V. However, it is 
evident from Figs. [4] and [5] that there are large matrix 
elements connecting such momenta to much larger mo- 
menta. This is directly associated with the repulsive 
core of the potential. For our discussion, we will adopt 
2 fm" 1 as the (arbitrary but reasonable) dividing line be- 
tween low and high momentum for nuclei. 

The consequences of the coupling to high momen- 
tum are readily seen in the probability density of the 
only two-body bound state, the deuteron. Consider the 
Argonne vis |5] curve in Fig. [6] The probability at 
small separations is significantly suppressed as a result 
of high-momentum components in the wave function. 
This suppression, called "short-range correlations" in 
this context, carries over to many-body wave functions 
and greatly complicates basis expansions. For exam- 
ple, in a harmonic oscillator basis, which is frequently 
the choice for self-bound nuclei because it readily al- 
lows removal of center-of-mass contamination, conver- 
gence is greatly slowed by the need to accommodate 
high-momentum components. The factorial growth of 
the basis size with the number of nucleons then greatly 
limits the reach of calculations. 
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Figure 6: Short-range correlations in the deuteron 3 Si probability den- 
sity from the AVI 8 potential (solid line). They are essentially elimi- 
nated by the RG evolution to lower flow parameter A (see Section [3T} . 



The underlying problem is that the resolution scale 
induced by the potential is mismatched with the basic 
scale of the low-energy nuclear states (given, for exam- 
ple, by the Fermi momentum). A solution is to elim- 
inate the coupling to high momentum. This is read- 
ily accomplished for spatial images (i.e., photographs) 
by Fourier transforming and then applying a low-pass 
filter — simply set the short wavelength parts to zero — 
and then transforming back. Let's try that for our 
Hamiltonian by setting to zero all of the matrix elements 
in Fig.|5]for k > 2 fm -1 . We test the implications by see- 
ing the effect on scattering phase shifts in the region of 
laboratory energies corresponding to k < 2fm _1 . 

The result is shown in Fig. [7] It would be unsurpris- 
ing that our filtered Hamiltonian fails close to the cut- 
off, but it is evident that there is a failure at all energies. 
What happened? The basic problem is that low mo- 
mentum and high momentum are coupled when solv- 
ing quantum mechanically for observables. For exam- 
ple, consider perturbation theory for the (tangent) of the 
phase shift (represented schematically here): 



<*|V|*> + J] 



(k\VW)(kf\V\k) 
(k 2 - k' 2 )/m 



+ •• 



(2) 



where a low momentum k is mixed with all other mo- 
menta at second order to a degree based on the size of 
the off-diagonal matrix elements. (As a computational 
aside, although momentum is continuous in principle, 
in practice we work on a discrete grid. This means that 
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Figure 7: Effect of a low-pass filter on observables: the ! So phase 
shifts. Note that the AVI 8 phase shifts reproduce experimentally ex- 
tracted phase shifts in this energy range. 



Eq. ^ becomes a matrix equation. For two-body po- 
tentials, roughly 100 x 100 matrices are sufficient.) The 
phase shift even at low energy or k will get large con- 
tributions from high k' if the coupling matrix elements 
{k\V\k') are large. 

How can we fix this? Our solution is to use a (short- 
distance) unitary transformation of the Hamiltonian ma- 
trix to decouple low and high energies. That is, insert 
the operator U^U = 1 repeatedly: 

E n = OWP = Wn\U*)UHU*UWn)) 

EE WnWVn) • (3) 

In doing so we have modified operators and wavefunc- 
tions but observables (measurable quantities) are un- 
changed. An appropriate choice of the unitary trans- 
formation can, in principle, achieve the desired decou- 
pling. This general approach has long been used in nu- 
clear structure physics and for other many-body applica- 
tions. The new feature here is the use of renormalization 
group flow equations to create the net unitary transfor- 
mation via a series of infinitesimal transformations. 

The renormalization group is well suited to this pur- 
pose. The common features of RG for critical phenom- 
ena and high-energy scattering are discussed by Steven 
Weinberg in an essay in Ref. |6|. He summarizes: 

"The method in its most general form can I 
think be understood as a way to arrange in 
various theories that the degrees of freedom 
that you're talking about are the relevant de- 
grees of freedom for the problem at hand." 



This is the essence of what is done with the low- 
momentum interaction approaches considered here: ar- 
range for the degrees of freedom for nuclear structure 
to be the relevant ones. This does not mean that other 
degrees of freedom cannot be used, but to again quote 
Weinberg [6 ]: "You can use any degrees of freedom you 
want, but if you use the wrong ones, you'll be sorry." 

The consequences of using RG for high-energy (par- 
ticle) physics include improving perturbation theory, 
e.g., in QCD. A mismatch of energy scales can gener- 
ate large logarithms that ruins perturbative convergence 
even when couplings by themselves are small. The RG 
shifts strength between loop integrals and coupling con- 
stants to reduce these logs. For critical phenomena in 
condensed matter systems, the RG reveals the nature 
of observed universal behavior by filtering out short- 
distance degrees of freedom. We will see both these 
aspects in our calculations of nuclear structure and re- 
actions. The end result can be said to make nuclear 
physics look more like quantum chemistry, opening the 
door to a wider variety of techniques (such as many- 
body perturbation theory) and simplifying calculations 
(e.g., by improving convergence of basis expansions). 

2.3. Summary points 

Low-energy nuclear physics is made difficult by a 
mismatch of energy scales inherent in conventional phe- 
nomenological potentials and those of nuclear structure 
and reactions. The renormalization group offers a way 
out by lowering the resolution, which means decoupling 
low- from high-momentum degrees of freedom. 

3. Overview of flow equations 

In Fig. [8] we show schematically two options for how 
the RG can be used to decouple a Hamiltonian matrix. 
The more conventional approach on the left in Fig. [8] 
lowers a cutoff A in momentum in small steps, with 
the matrix elements adjusted by requiring some quan- 
tity such as the on-shell T-matrix to be invariant. (In 
practice this may be carried out by enforcing that the 
half-on- shell T-matrix is independent of A.) Matrix el- 
ements above A are zero and are therefore trivially de- 
coupled. When adapted to low-energy nuclear physics, 
this approach is typically referred to as "Vi ow ^" EEL 

The approach we will focus on is illustrated in Fig. [8] 
on the right, in which the matrix is driven toward band- 
diagonal form, achieving decoupling again but without 
truncating the matrix. The corresponding RG was de- 
veloped in the early 1990's by Wegner (HlEHTOl for con- 
densed matter applications under the name "flow equa- 
tions" and independently by Glazek and Wilson ifTTll for 
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Figure 8: Schematic Vi 0W fc RG evolution (left) and flow equation RG evolution (right). 
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Figure 9: SRG flow of the AVI 8 NN potential in the 3 Si channel at selected values of the flow parameter A 0. 



solving quantum chromodynamics in light- front formal- 
ism under the name "similarity renormalization group" 
(SRG). Only in the last five years was it realized that 
the band-diagonal approach is particular well suited for 
low-energy nuclear physics, where it is technically sim- 
pler and more versatile than other methods [Q21 |2) . We 
will apply formalism closer to the flow equation formu- 
lation, but generally use the simple abbreviation SRG. 
We have introduced a cutoff-like parameter A that serves 
as a momentum decoupling scale. Elsewhere we will 
also use the natural flow parameter s = 1 /A 4 . 



3.1. Flow equation basics 

The basic flow equation is a set of coupled differen- 
tial equations for the (discrete) matrix elements of the 
Hamiltonian matrix (or the potential in practice, because 
we fix the kinetic energy matrix by construction). Let's 
see an example in action before stepping back and con- 
sidering details. For the two-body potential in a partial- 



wave momentum basis, the flow equation takes the form 

^(k,k')oc-(e k -e k ,) 2 V A (k,k') 
dX 

+ Y}& + ^ - ItqWAik, q)V A (q, k') , (4) 

q 

where k and k' are the relative momenta and e^ = 
H 2 k 2 /M and we have omitted an inessential constant. 
The evolution is continuous, but snapshots at selected 
A values are shown in Fig. [9] (Note that the axes are the 
kinetic energy k 2 .) The evolution toward diagonaliza- 
tion is evident, with the width of the band in k 2 roughly 
given by A 2 . 

It is evident that the off-diagonal matrix elements are 
driven toward zero, so we expect that a low pass filter 
will now be effective. Indeed it is, as shown in Fig. 10 
where NN phase shifts for the AVI 8 potential in a va- 
riety of channels are compared with the results from 
applying a low-pass filter to the original (dotted) and 
evolved (dashed) potentials. The evolved result agrees 
up to the low-pass cutoff. Note that if this cutoff is not 
applied, the phase shifts for the evolved Hamiltonian 
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Figure 10: Nucleon-nucleon (NN) phase shifts showing the effect of a low-pass filter at & max = 2.2fm _1 in various partial-wave channels I13II14I . 
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The dashed line is the results when the potential is evolved first and then cut at k max . 



agree precisely at all energies, because the transforma- 
tion is unitary in the two-body system (and this is pre- 
served to high numerical precision). 

If we now revisit the consequence of a repulsive core 
for the deuteron probability density, we see in Fig. [6] 
that the "wound" in the wave function at small r is filled 
in as the core is transformed away. (Note: it may look 
like the normalization is not conserved, but if we multi- 
plied by r 2 , we would see that the area under the curves 
is the same, as are the large-r tails.) Thus the short- 
range correlations in the wave function are drastically 
reduced. This means that the physics interpretation of 
various phenomena is altered, even though the observ- 
able quantities such as energies and cross sections are 
unchanged. (The long-range part of the wave function 
is preserved; this is related to the asymptotic normal- 
ization constants, which can be extracted from experi- 
ment.) We cannot immediately visualize the changes in 
the potential in coordinate space in a conventional plot 
like Fig. [3] however, because it is now non-local, which 
is to say it is not diagonal in coordinate representation: 



V(T)l/f(T) 



/■ 



d 3 r'V(r,r')<Kr') 



(5) 



This is a technical problem for certain quantum many- 
body methods (such as Green's function Monte Carlo) 
but not for methods using harmonic oscillator matrix el- 



ements. 

We can visualize the evolution approximately by con- 
sidering a local projection of the potential: 



V A (r) 



■/■ 



d 3 r'V A (r 9 r') 



(6) 



which leaves a local V(r) unchanged. For a non-local 
potential, this roughly gives the action of the potential 
on long- wavelength nucleons. This is shown in Figs. 11 



and 12 where in the former we see the central potential 
dissolving and in the latter similar effects on the ten- 
sor part of the potential 1 15 ]. Also evident is the flow of 
the two potentials, initially quite different (potentials are 
not observables!), toward a universal flow at the lower 
values of X. Very recent work suggests that such a local 
projection may capture most of the physics of the full 
evolved potential, with the effects of the residual poten- 
tial calculable in perturbation theory. This may open 
the door to using RG-evolved potential with quantum 
Monte Carlo methods. 

The consequences of low-momentum potentials for 
harmonic-oscillator-basis calculations are illustrated in 



Figs. 13 and 14 (We'll refer to the method used, which 
is a direct diagonalization of the Hamiltonian matrix and 
therefore is variational, as no-core full configuration, or 
NCFC.) The original potential in this case is already 
soft (that is, there is much less coupling to high mo- 
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mentum than in the AVI 8 potential), but convergence 
in a harmonic oscillator basis with A^ max shells for ex- 
citation is slow (these are the "Original" curves). Note 
that the matrix dimension grows rapidly with N max and 
the number of nucleons A; for example A^ max = 8 has 
dimension about 50,000 for 4 He but over one million 



for b Li (see Fig. 15 for other examples). But with SRG 



evolution, there is vastly improved convergence. The 
convergence is also smooth, which makes it possible 
to reliably extrapolate partly converged results to the 
A/max — > oo limit. Nevertheless, the rapid growth of 
the basis with A is still a major hindrance to calculat- 
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Figure 15: Dimension of the Hamiltonian matrix for NCFC/NCSM 
calculations as a function of the harmonic oscillator basis size (iVmax 
shells) for selected nuclei. 



ing larger nuclei. One solution being explored by Roth 
and collaborators (HDDS) is to use importance sampling 
of matrix elements, in which only a fraction of the full 
matrix is used. This technique is enabled by the RG 
softening of the potential, which allows the importance 
to be evaluated perturbatively. 

Let's now return to the basics of SRG flow equations. 
We wish to transform an initial hamiltonian, H = T + V 
in a series of steps, labeled by the flow parameter s: 



H s = U S HU] = T + V S 



with 



UlUs = u s u 



1 



(7) 



(8) 



Note that the kinetic energy T is taken to be independent 
of s. Differentiating with respect to s\ 



dH s dU s + + 4. 

s - - -77' 77 HIP + I J HIP J J - 



= [71s, H s ] 



with 



dV, rf 



(9) 



(10) 



The anti-Hermitian generator rj s can be specified by a 
commutator of H s with a Hermitian operator G s : 

T] S = [G S ,H S ], (11) 

which yields the flow equation (with T held fixed!), 



dH s dV s 
ds ds 



(12) 



The operator G s determines the flow and there are many 
choices one can consider. 

Probably the simplest example we can consider is just 
a two-state system |20|. Let H = T + V, where 



T\i) = £i \i) and Vij = (i\V\j) 
Then we can choose G v = T and 



dHs 

ds 



[[T,H S ],H S ] 
becomes (with the s dependence implicit) 

-(S-e/Vy 

+ J^iei + ej - 2e k )V ik V kj . 



(13) 



(14) 



ds V,J 



(15) 



For a two-level system with i = {a, b), we can express 
the flowing Hamiltonian in terms of Pauli matrices: 



T = -(e a + e b )l + -(e a - e b )cr z 



(16) 



and 



V s = ^{Vaa + Vbb)I 



+ j(V. 



V bb )cr z + V ab cr x 



(17) 



The solution to the SRG flow equation is easily found. 
It is convenient to parametrize the result in terms of 6(s) 
with constant qj: 



d6 
ds 



-2(e a - e b )aj sin 0(s) 



(18) 



<D 4 



1 1 1 1 1 1 1 1 
8 =3, 8 =9 


^ 


yS -< aa 


^^ V K 


/ 2-body matrix elements 


i i i i i i i i 



3 

A,= l/s 



4 

1/4 



10 



Figure 16: SRG flow for a simple two-state system (see text). 
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Figure 17: SRG flow in stages. On the far left and far right are potentials in the l So channel evolved to A = 2.0 fm ! and A = 1 .5 fm l , respectively. 
The middle panels show the first (left) and second (right) terms on the right side of Eq. (23). 



where 



0(s) = 2 tan" 1 [t<m(6(0)/2)e- 2(6a - €b)cos ] , (19) 
co cos 6 = (e a - e b + V aa - V bb )/2 , (20) 



and 



cosinO = V a b • (21) 

The resulting flow is plotted for sample energies e a and 



€b in Fig. 16 We clearly see the off-diagonal matrix ele- 
ment V a b driven to zero. (Try reproducing this in Math- 
ematical) 

For a nuclear two-body (NN) potential in a partial- 
wave momentum basis with tj s = [T, H s ], we project on 
relative momentum states \k) using 1 = | C° \q)q 2 dq{q\ 
with H 2 /M = 1. The flow equation reduces to: 



dVs_ 
ds 



= [[Trt,V s ],H s ] with T rQl \k) = e k \k) (22) 



and A 2 = 1/ Vs. T YQ \ is the relative kinetic energy of the 
nucleons. Then 

^(*,*0 oc -(6 k -6 k ,fV X (Kk') 

dA 

+ Yj& + ^ - teqWAik, q)V A (q, k') . (23) 



This particular equation is for A = 2, but the results are 
generic if one lets k represent a set of Jacobi momenta. 
The first term in Eq. ([23]) drives Vx toward the diagonal: 



V A (k,k?) = V A=00 (k,V)e 



-[{e k -e k >)IA 2 



(24) 



which can be visualized in Fig. 17 These panels repre- 
sent a sequence from A = 2.0 to A = 1.5. The potentials 
at the beginning and the end are on the outside, while 
the two middle panels are the first and second term of 
Eq. ( [23] ). For off-diagonal matrix elements, the first term 
is numerically dominant and each element is driven to 
zero according to Eq. ( [24] ), with further off-diagonal ele- 
ments changing more rapidly. Note that the width of the 
diagonal is given roughly by A 2 , in accord with Eq. p4|. 



A more general proof follows if we use the genera- 
tor advocated by Wegner, which includes the diagonal 
part of the Hamiltonian, H d . Call the diagonal elements 
H it = e i9 then (with rj s = [H d , H s ]), 



dHtj 
ds 



But 



= (i\[[H d ,H 5 ],H s ]\j) 

= ^j(ei + ej-2e k )H ik H kj 

k 

= 2Y J (e i -e k )\H ik \ 2 . 



*E^ = 2 Z^ 



(25) 



ds 



i 



= 4 J] «.-(«.- - e k )\H ik \ 2 

i*k 

= 2Y J (ei-e k ) 2 \H ik \ 2 >0. (26) 

frk 

Now use: 

TrH 2 = const. = J] IH^ 2 = J] IH/ + J] |Hij| 2 , 

(27) 



y 



i*j 



to obtain 



ds 



&j 



ds 

i 

i+k 



e, - e k f\H ik \ 2 < . (28) 



Thus, in the absence of degeneracies, the off-diagonal 
matrix elements will decrease (or at least remain un- 
changed) 1 1 ]. 

This feature of the diagonal generator is desirable, but 
we also note that for nuclear Hamiltonians in a momen- 
tum basis, the diagonal is completely dominated by the 
kinetic energy, so Hj « T YQ \ is a very good approxi- 
mation. Can this break down? Glazek and Perry lf2TI 
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showed that it can (see also Wendt et al. l22l ). Re- 
consider the proof of diagonalization, but now with 
G s = T rQ \. Now we have Ha = et and T m \\i) = 6i\i), 
so that 

dHu 
ds 

So consider 



= 2Y J (E i -e k )\H ik \ 2 . 



(29) 



d X l^l 2 = 2^^^ = 4 Z^ " *>l»ftl : 



ds 



/V/v 



(30) 



from which we conclude 



^ J) |ify| 2 = -2 2^ - e k )(€t - e k )\Hik\ 2 • (31) 



ds 



i*3 



itk 



Thus the off-diagonal decrease depends on e t - e k « 
6/ - £&. But there is the possibility of this not being true, 
e.g., if there are spurious bound states as in large-cutoff 
EFT [221. 



3.2. Alternative generators 

Other choices of generator are also possible. Recent 
work by Shirley Li while an undergraduate physics ma- 
jor at Ohio State explored choices designed to accelerate 
evolution. In particular, one can choose G s as 



G s = - 



A 1 



1 + TWA 2 



or 



Gs = _aV^ /a2 



c + r rd + • • 



• c + r rd + • 



(32) 



(33) 



The expansions show that these reduce to the conven- 
tional r re i for momenta that are small compared to the 
cutoff parameter A (not to be confused with A). For 
A = 2 fm" 1 , the low energy part of the potential is still 
decoupled but there is much less evolution at high en- 
ergy, which makes it computationally much faster and 
allows evolution to low A. See Ref. l23l for details. 

The usual approach to Vi ow fc RG evolution (left dia- 
gram in Fig. [8]) is based on the Lippmann-Sch winger 
equation for the half-on- shell T-matrix illustrated in 
Fig. 



18 



A cutoff A is imposed on the integral in the 
second term and we demand that matrix elements of T 
are invariant with an infinitesimal reduction of A. That 
is, we require dT(k,k';Ek)/dA = 0, which establishes 
an RG equation for V\- In contrast to the SRG equation, 
which is second order in the running potential, the Vi ow & 





Figure 18: Schematic version of the Lippmann-Sch winger equation 
for the T-matrix, with cutoff A on the intermediate states. The Vi ow fc 
potential Va is determined by requiring half-on- shell matrix elements 
of this equation to be invariant under changes in A (7] . 



RG equation has the T-matrix on the right side. Thus 
T(k\k;k 2 ) = V NN (k\k) 

p 2 dp 



2_ C 
-9 

n Jo 



k 2 -p 2 



^V m (k\p)T(p,k;k 2 ) 

= VL*(*\*) 
_2 p r^Vt k (k\ P )T( P ,k;k 2 )^ 2 

n Jo 



k 2 -p 2 



p dp 
(34) 



for all k, k! < A. (Note: we are using standing- wave 
boundary conditions for numerical reasons; this is often 
called the K-matrix.) From dT/dA = 0, we get the 
Viow^RG equation: 



^vL k (^ = 



2Vt,(k\A)T A (A,k;A 2 ) 



1 - (k/A) 2 



(35) 



Note that the full T matrix appears on the right side, 
in contrast to the partial- wave SRG flow equation (with 

G s - T re i), 

^V A (k 9 kf) cc -(E k -e k ,) 2 V A (k,k') 
dX 



+ 2>* + 6 *' " 26 «) W 9)VM k') , (36) 



which is only second-order in the potential. We can also 
define smooth regulators for V\ ow k as in Ref. l24l . 

The evolution of NN potentials using the Vi ow fc 
method is illustrated for a sharp and smooth regulator in 
Fig. [19] Comparing the Vi ow & flow in these figures to the 



SRG flow in Fig. [9] we see the same decoupling of low 
and high momentum. Other non-RG unitary transfor- 
mations (which perform the transformation all at once, 
rather than in steps) also decouple; an example is the 
UCOM method |25 ], which has close connections to the 
SRG (see Refs. [26,27]). 
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Figure 19: V\ ov/ k flow of AVI 8 in the 3 Si channel for a sharp (top) and smooth (bottom) regulator |24|. 
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Figure 20: SRG sharp block-diagonal flow of AV18 in the l P\ channel at A = 4, 3, 2, and 1 fm -1 [28]. The initial N 3 LO potential is from Ref. (29] 
and A = 2 fm -1 . The axes are in units of k 2 from to 1 1 fm -2 and the color scale is from -0.5 to +0.5 fm. 



It is also possible to choose a generator that repro- 
duces the block diagonal (as opposed to band diagonal) 
form of the V\ ow k RG shown schematically in Fig. [8] 
except that the transformation will be unitary. In partic- 
ular, we can use 



^ = [[G S9 H S IH S ] 
as 



with 



PH S P 






QH S Q 



(37) 



(38) 



where projection operators P and Q = 1 - P are simply 
step functions at a given A in partial-wave momentum 
representation. An example of the subsequent flow is 
shown in Fig. 20 for the ^i channel |28|. To get the 



fully block-diagonal form, one would have to evolve to 
A - oo. But in practice, going to A = 1 fm -1 is sufficient 
for essentially complete decoupling at A = 2 fm -1 . 



The proof of block diagonalization (see Gubankova et 
al. (30l[3Tl) goes as follow. The generator r] s = [G s , H s ] 
is non-zero only where G s is zero, which means in the 
off-diagonal blocks. This will then evolve the potential 
in this same pattern (this is generically true, if one de- 
sires a different pattern [28]). A measure of off-diagonal 
coupling of H s is 

Tr[(QH s P)\QH s P)] = Tr[PH s QH s P] > . (39) 

Now we can calculate the derivative of this expression 
by applying the SRG equation ( [37] ): 



ds 



Tv[PH s QH s P] = 



Tr[Pr, s Q(QH s QH s P - QH S PH S P)] 
+ Tr[(PH s PH s Q - PH s QH s Q)Qr] s P] 
= -2Tr[(Qr 1s P)\Qr 1s P)]^0. (40) 
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Thus the off-diagonal QH S P block will decrease as s 
increases. 

The low-momentum block of this SRG is found to be 
remarkably similar to the corresponding Vi 0W £ RG po- 
tential. Examples are shown in Figs. [21] and [22] How- 
ever, proving that these different RG approaches should 
yield the same potential remains an open problem. 



(a) k 1 (fm" 1 ) 




n0.5 



(fm) 



I 



0.5 



Figure 21: Comparison of momentum- space (a) V\ ow k and (b) SRG 
block-diagonal 3 Si potentials with A = 2fm _1 [28]. 
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Figure 23: Convergence of three-body energies for two potentials and 
several sets of A's. 
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Figure 22: Same as Fig. [21] only a surface plot. 
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3.3. Many-body forces 



In Fig. 23 the convergence of the triton ground- state 
energy with harmonic oscillator basis size (N mdLX Haj ex- 
citations) is shown for two chiral EFT NN potentials 
and for the corresponding SRG-evolved potentials at 
/i's from 4 far 1 to lfm -1 . In accord with our previ- 
ous discussion, we see increasingly rapid convergence 
as A decreases. (Note that there is softening already 
at A = 3 far 1 for the N 3 LO EFT with A = 600 Me V, 
which corresponds to 3 fm -1 . The moral is that it is not 
sufficient to simply compare the cutoff numerically to 
the decoupling scale A.) However, we also see that the 
converged (or extrapolated) binding energies are differ- 
ent for each A\ How could this be, if the SRG is sup- 
posed to generate unitary transformations? 

There are further signs of trouble, such as nuclear 
matter failing to saturate after the two-body potential is 
softened by RG evolution. Let's review the facts about 



energy per particle for pure neutron matter (Z = 0) and 
symmetric nuclear matter (N = Z = A/ 2) from cal- 
culations adjusted to be consistent with extrapolations 
from nuclei. Neutron matter has positive pressure, while 
symmetric nuclear matter saturates; that is, there is a 
minimum at a density of about 0.16 fm 3 with a binding 
energy of about -16MeV/A. Reproducing this mini- 
mum with microscopic interactions fit only to few-body 
data is one of the holy grails of nuclear structure the- 
ory. But if we evolve an NN potential with either the 
Viow£ RG or the SRG, we find that nuclear matter does 
not converge. This is shown for Vi 0W £ by the "NN only" 



curves in Fig. 25 and similar behavior is found for SRG. 
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Figure 24: Binding energies of nuclear and neutron matter from Ak- 
mal et al. for several equations of state |32|. 



uniform nuclear matter. Figure 24 shows the binding 
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Figure 25: Nuclear matter at third order in perturbation theory in the 
particle-particle channel (this appears sufficient for convergence) us- 
ing a V\ 0W k RG-evolved potential |33]. The "NN only" curves include 
no three-body interactions while the "NN + 3N" curves include a 3NF 
fit to the triton binding energy and the alpha particle radius. 



The failure of softened low-momentum potentials to 
reproduce nuclear matter saturation should sound famil- 
iar to anyone who knows the long history of low-energy 
nuclear theory. There were active attempts to transform 
away hard cores and soften the tensor interaction in the 
late sixties and early seventies. But the requiem for soft 
potentials was given by Bethe in 1971 l34l : "Very soft 
potentials must be excluded because they do not give 
saturation; they give too much binding and too high den- 
sity. In particular, a substantial tensor force is required." 
The next thirty-five years were spent struggling to solve 
accurately with such "hard" potentials. But the story is 
not complete: the three-nucleon forces (3NF) were not 
properly accounted for! 

Three-body forces between protons and neutrons 
have a classical analog in tidal forces: the gravitational 
force on the Earth is not just the pairwise sum of Earth- 
Moon and Earth-Sun forces. Quantum mechanically, an 
analog is with the three-body force between atoms and 
molecules, which is called the Axilrod-Teller term and 
dates from 1943 ll35l . The origin is from triple-dipole 
mutual polarization. It is a third-order perturbation cor- 
rection, so the weakness of the fine structure constant 
means that these forces are usually negligible in metals 
and semiconductors. However, in solids bound by van 
der Waals potentials it can be significant; for example, it 
is ten percent of the binding energy in solid xenon [36|. 
This is the same relative size as needed in the triton. 

In general, three-body forces arise from eliminating 
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Figure 26: Sources of internucleon three-body forces. On the left, the 
intermediate state includes an excitation of the nucleon (A or N*). 
On the right, the intermediate state includes a virtual anti-nucleon 
(N). When these degrees of freedom are integrated out, the result- 
ing nucleons-only vertex is a 3NF. 



degrees of freedom. In the nuclear case, this can mean 
eliminating excited states of the nucleon (N* or A) or 
from relativistic effects; see Fig. 26 for diagrammatic 
representations. If the intermediate states are not in- 



cluded in the low-energy degrees of freedom, we have 
irreducible vertices with three nucleons in and three nu- 
cleons out. But 3NF's also result from the decoupling 
of high-momentum intermediate states, whether they 
are eliminated explicitly by a cutoff (as with the Vi ow & 
RG) or the degree of coupling modified (as with the 
SRG). Omitting three-body forces leads to model de- 
pendence: observables will depend on the decoupling 
scale, whether it is the Vi 0W £ A or the SRG X. This de- 
pendence also becomes a tool, because it is a diagnostic 
for errors (more on this later). To eliminate this depen- 
dence, the 3NF at different A or X must be either fit or 
evolved. 

It is easy to see that RG flow equations lead to many- 
body operators. Consider the SRG operator flow equa- 
tion written with second-quantized a's and o^'s: 



dVs_ 

ds 



\\/ i CL^a , y k a^a^aa |, y k cficJaa \ 
G * 2-body 2-body 

y cfi a* a* aaq + • • • (41) 

3-body! 



= ... + 



where the second equality reflects that even if the initial 
Hamiltonian is two-body, the commutators give rise to 
three-body terms. (For future reference, recall that the 
creation and destruction operators are always defined 
with respect to a single-particle basis and a reference 
state, which in this case is the vacuum.) 

As the evolution continues, there inevitably will be 
A-body forces (and operators) generated. Is this a prob- 
lem? Not if these "induced" many-body forces are the 
same size as those that naturally occur. Indeed, nuclear 
three-body forces are already needed in almost all po- 
tentials in common use to get even the triton binding 



R.J. Furnstahl / Nuclear Physics B Proceedings Supplement 00 (2013) l-\37\ 



15 





Ci,C 3 ,C 4 

Figure 27: Leading three-body forces from chiral EFT. These contri- 
butions represent three different ranges: long-range 2-pion exchange, 
short-range contact with one-pion exchange, and pure contact interac- 
tion. 



energy correct. In fact, low-energy effective theory tell 
us generalized diagrams such as those in Fig. [26] with 
four or more legs imply that there are A-body forces 
(and operators) initially! 

However, there is a natural hierarchy predicted from 
chiral EFT, whose leading contributions are given in 



Fig. [27] (we'll return to this in Section |47T| and supply 
additional details). If we stop the flow equations be- 
fore induced A-body forces are unnaturally large or if 
we can tailor the SRG G s to suppress their growth, we 
will be ok. (Another option is to choose a non-vacuum 
reference state, which is what is done with in-medium 
SRG, to be discussed later.) Note that analytic bounds 
on A-body growth have not been derived, so we need to 
explicitly monitor the contribution in different systems. 
But the bottom line that makes the SRG attractive as 
a method to soften nuclear Hamlitonians is that it is a 
tractable method to evolve many-body operators. 

To include the 3NF using SRG with normal-ordering 
in the vacuum, we start with the SRG flow equation 
dHJds = [[G S ,H S ],H S ] (e.g., with G s = T rel ). The 
right side is evaluated without solving bound-state or 
scattering equations, unlike the situation with Vi 0W yb so 
the SRG can be applied directly in the three-particle 
space. The key observation is that for normal-ordering 
in the vacuum, A-body operators are completely fixed in 
the A-particle subspace. Thus we can first solve for the 
evolution of the two-body potential in the A = 2 space, 
with no mention of the 3NF (either initial or induced), 
and then use this NN potential in the equations applied 
to A = 3. 

What about spectator nucleons? There is a decou- 
pling of the 3NF part. We can see this from the first- 
quantized version of the SRG flow equation, 

ds 



dV i2 dVj 3 dV 2 3 dVn3 



ds 



ds 



ds 



ds 



[[T mh V s lH s ]., 



(42) 



where we isolate the contributions from each pair and 
the 3NF. Using each SRG equation for the two-body 
derivatives, we can cancel them against terms on the 
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Figure 28: Triton binding energy during SRG evolution. The three 
curves are for an initial potential with only NN components where 
the induced 3NF is not kept ("NN-only"), for the same initial NN 
potential but keeping the induced 3NF ("NN + NNN-induced"), and 
with an initial NNN included as well ("NN + NNN"). 



right side. The result is 1 12]: 

dV U 3 



ds 



= [[ri2,Vl2],(r 3 + Vl3 + V23 + Vi23)] 

+ {123-> 132} + {123 -^ 231} 
+[[T Tch V 12 3lH s ]. (43) 

The key is that there are no "multi- valued" two-body 
interactions remaining (i.e., dependence on the excita- 
tion energy of unlinked spectators); all the terms are 
connected. An implementation of these equations in a 
momentum basis would be very useful and has very re- 
cently been achieved by Hebeler [37]. But an alternative 
approach has also succeeded: a direct solution in a dis- 
crete basis (38j EH HZ). 

The idea is that the SRG flow equation is an opera- 
tor equation, and thus we can choose to evolve in any 
basis. If one chooses a discrete basis, than a separate 
evolution of the three-body part is not needed. This was 
first done for nuclei by Jurgenson and collaborators in 
2009 using an anti- symmetrized Jacobi harmonic oscil- 
lator (HO) basis Q6). The technology for working with 
such a basis had already been well established for appli- 
cations to the no-core shell model (NCSM) (39). This 
approach leads to SRG-evolved matrix elements of the 
potential directly in the HO basis, which is just what 
is needed for many -body applications such as NCFC or 
coupled cluster. 

the comparison of two-body-only to full 



In Fig. 28 



two-plus-three-body evolution is shown for the triton 
( 3 H). The NN-only curve uses the evolved two-body 
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Figure 29: Alpha particle binding energy during SRG evolution. The 
curves correspond to those in Fig. [28] 



Figure 3 1 : Contributions to the triton binding energy during SRG evo- 
lution. Plotted are the expectation values of the kinetic energy, the 
two-body potential, and the three-body potential 1171 . 
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Figure 30: Correlation plot of the binding energies of the alpha parti- 
cle and the triton. The dotted line connects (approximately) the locus 
of points found for phenomenological potentials, and is known as the 
Tjon line (40). 



potential. The change in energy with X reflects the vi- 
olation of unitarity by omission of the induced three- 
body force. When this induced 3NF is included ("NN 
+ NNN-induced"), the energy is independent of X for 
A = 3. If we now turn to the alpha particle ( 4 He) in 



Fig. 29 we see similar behavior, except now the inclu- 
sion of the induced 3NF does not lead to a completely 
flat curve at the lowest X values. If there is sufficient 
convergence, this is a signal of missing induced 4NF. 

In both cases, it is evident that starting with an initial 
NN-only interaction (in this case, an N 3 LO(500MeV) 



third line in each plot of Figs. [28] and [29] shows that an 
initial 3NF (labeled NNN) contribution leads to a good 
reproduction of experiment. The triton energy is part of 
the fit of this initial force, but the alpha particle energy 
is a prediction. Note that the magnitude of the NN-only 
variation is comparable to the initial 3NF needed. This 
is an example of the natural size of the 3NF being mani- 
fested by the running of the potential (which is, in effect, 
the beta function). 

The nature of the evolution is illustrated in Fig. 30 
which is a correlation plot of the binding energies in 
each nucleus. The dotted line is known as the Tjon 
line for NN-only phenomenological potentials. It was 
found that different potentials that fit NN scattering data 
gave different binding energies, but that they clustered 
around this line. With the SRG evolution starting with 
just an NN potential, the path follows the line, passing 
fairly close to the experimental point. With an initial 
NNN force and keeping the induced 3 -body part, the 
trajectory is greatly reduced (see inset), at least until X 
is small. 



Figures [31] and [32] show individual contributions to 
the energy in the form of ground- state matrix elements 
of the kinetic energy, two-body, three-body, and (im- 
plied) four-body potentials. The hierarchy of contri- 
butions is quite clear but the graphs also manifest the 
strong cancellations between the NN and kinetic energy 
contributions. These cancellations magnify the impact 
of higher-body forces. Even so, it appears that a trunca- 
tion including the NNN but omitting higher-body forces 
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is workable, particularly with A > 1.5 fm -1 . But what 
about the A dependence of the 4NF (and beyond)? 
This A dependence is the topic of current research. In 
results for 6 Li are shown [ 17 ] 



Fig. 33 results for °Li are shown [17]. Assessing these 
results is made difficult because of insufficient conver- 
gence of comparison calculations (the shaded areas) and 
of calculations at the large A values. Note, however, that 
the variations are with the 1 MeV level; nevertheless we 
expect to do better. 

Roth and collaborators have since used importance 
truncated NCSM (IT-NCSM) to extend these calcula- 
tions to much higher N max and all the way to 16 1 19]. 
They find that the SRG with initial NN-only but includ- 
ing the induced NNN shows only small running with 
A. On the other hand, with increasing A they find sig- 
nificant deviations. This has been traced to the influ- 
ence of the initial long-range NNN interaction. If they 
lower the cutoff of this part of the interaction, then ap- 
proximate SRG unitarity is restored and with coupled- 
cluster methods they find reasonable results even for 
medium-size nuclei (although not yet with fully consis- 
tent Hamiltonians) flTTl . 

3.4. Summary points 

Renormalization group flow equations dramatically 
reduce correlation in many-body wave functions, lead- 
ing to faster convergence of many-body calculations. 
Flow equations (SRG) achieve this lower resolution 
by decoupling via a series of unitary transformations, 
which leave observables invariant (if no approximations 
are made) but alter the physics interpretation. Few-body 



forces are inevitable, but the flow-equation approach al- 
lows the evolution of vacuum interactions. 



4. Features of SRG applied to nuclear problems 

4.1. Chiral EFT, many -body forces, and the SRG 

Before continuing with the SRG for nuclear systems, 
let's say a bit more about chiral effective field theory 
(EFT). In the SRG flow equations, the input interaction 
is merely an initial condition; the equations are the same 
whether we start with an EFT potential or a more phe- 
nomenological potential such as Argonne vig. However, 
increasingly nuclear theorists are moving toward using 
EFT interactions because they promise a more system- 
atic construction of many-body forces and consistent 
operators. 

There are three fundamental ingredients of an effec- 
tive field theory (e.g., see Ref. |42|). The first is to use 
the most general lagrangian with low-energy degrees of 
freedom consistent with the global and local symme- 
tries of the underlying theory. For nuclei, the underly- 
ing theory is quantum chromodynamics (QCD). We can 
identify a hierarchy of nuclear QCD scales: 



Mqcd ~ 1 GeV [M hadrons , AnfA 



• M m 



100M6V [k F ,f n ,m n ,6w] 



• ^nuc/^QCD ~ 10 MeV [nuclear binding energy 
per nucleon] 
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For all but the lightest nuclei, the EFT of choice at 
present draws a line between the first two levels to de- 
fine the high-energy and low-energy scales. This is the 
chiral EFT, with degrees of freedom consisting of the 
nucleon (proton and neutron) and the pion (43). In the 
near future, the A resonance will be included as well 
because of the small mass difference with the nucleon, 
San- Besides the usual space-time symmetries, terms 
in the chiral EFT lagrangian are constrained by the re- 
quirements of spontaneously broken (as well as explic- 
itly broken) chiral symmetry (43). 

The other two ingredients are the declaration of a reg- 
ularization and renormalization scheme, and the identi- 
fication of a well-defined power counting based on well- 
defined expansion parameters. The separation of scales 
provides the expansion parameter as ratios of Q/Mqqd, 
where Q is one of the quantities lumped together above 
as M nuc . The chiral EFT potentials used here are derived 
using a momentum cutoff and what is called "Wein- 
berg counting", in which the counting is done at the 
level of an irreducible potential that is summed non- 
perturbatively with the Lippmann-Sch winger equation. 
This scheme has been criticized because it does not al- 
low systematic renormalization (meaning, in this con- 
text, the removal of cutoff dependence at each order), 
which limits the range of cutoffs used. This in turn hin- 
ders the validation of the EFT, because the sensitivity 
to the cutoff is used as a measure of uncertainties, See 
Ref. 1 43 ] for a thorough overview of chiral EFT for nu- 
clei and the status of the power counting and renormal- 
ization controversies. 

With any scheme, however, the power counting im- 
plies a hierarchy of many -body contributions. Wein- 
berg counting associates a power v of Q/Mq C d with di- 
agrams for the potential, where 

y = -4 + IN + 2L + ^(d,- + m - 2) . (44) 

The definitions and details of its implementation can 
be found in Ref. |43) . For our purposes, the relevant 
term is "IN", which says that adding a nucleon to go 
from an A-body potential to an A + 1-body potential 
generally suppresses the contribution by Q 2 /Mqqd- In 
the theory without A's, the suppression of the leading 
3NF compared to the leading NN interaction is actually 
(2/Mqcd) 3 , and a four-body force first appears at or- 
der (2/Mqcd) 4 [43]. It is this hierarchy that we want to 
preserve as we run our SRG flow equations. 

As noted, the flow equation technology discussed 
here does not rely on a particular implementation of 
chiral EFT, except that the SRG is inherently non- 
perturbative. This apparently excludes alternative renor- 



malization and power counting schemes that require a 
perturbative treatment beyond leading order (4U |45l 
[46). The current belief is that the two approaches should 
give comparable results as long as the EFT cutoff is 
taken to be of order Mqcd, but the issue is far from set- 
tled 1431 . 

Another consideration is whether one could bypass 
the SRG by simply applying chiral EFT with a lower 
cutoff. Indeed, there exists low-cutoff implementations 
on the market that display similar characteristics to 
low-momentum RG-evolved interactions. However, the 
lower cutoff also means the effective expansion param- 
eter is smaller, and therefore the truncation error is re- 
duced. With the RG, one preserves the truncation error 
from the cutoff of order Mq C d- An additional advantage 
of the RG is the controlled variation of the decoupling 
scale, which provides a tool for assessing errors from 
the Hamiltonian truncation and many-body approxima- 
tions. But this is also not a settled issue and further 
study would be welcome. 

4.2. More perturbative nuclear systems flow equations 
Earlier we mentioned the role of RG in high-energy 
physics in improving perturbation theory. Much of low- 
energy nuclear physics is intrinsically non-perturbative 
because of large scattering lengths and bound states, so 
how do we quantify "perturbativeness"? We study the 
convergence of the Born series for scattering to see how 
this can be done. 

Consider whether the Born series for the T-matrix op- 
erator at a given (complex) z, 

T( Z ) = V+V V+V V V+... (45) 



H z-H z-H 

converges. This is something like a geometric series, 
for which we know clearly the convergence criterion: 
1 + w + w 2 + • • • diverges if \w\ > 1. We get a clue for 
how to use this if we consider a bound state \b) and the 
special value z = Eb, which is the bound- state energy. 
Then we can rearrange the Schrodinger equation 



(H + V)\b) = E b \b) 



to the form 



1 



-V\b) = \b) 



(46) 



(47) 



Eb - Hq 

and look at T(E b )\b). Using Eq. ( |47| ) repeatedly, the 
divergence is manifest, i.e., we get V(l + 1 + 1 + • • • )\b). 
Now we see that we can generalize Eq. ( |47| ) for fixed 
E by looking for eigenstates of (E - Hq)~ 1 V with eigen- 
value T]y(E), 



E-H ( 



■vir v > = i7 v (E)|r v >. 



(48) 



0"2 
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>n 



^ J00 



.> 



R.J. Furnstahl / Nuclear Physics B Proceedings Supplement 00 (2013) l-\37\ 
J \2 5 T2 



-3 -2 -1 
I IP h# — I 1 \-^+- 



19 



500 



300 



VAJ8 /^ 
HS!qd3 If 






- C0L6 / 1 
"L6bn|2!Ael J 




1 

it 





,3° 


C[J9IJU6| ■ 

1 1 1 1 1 1 



Figure 34: The potentials of Fig. [3] inverted as part of the Weinberg 
eigenvalue analysis (see text). 



Then with T applied to these eigenstates, there is mani- 
festly a divergence for \rj v (E)\ > 1: 



r(£)|r v > = v|r v >(i + 77 v + rf y + •••)• 



(49) 



So we characterize the perturbativeness of a potential 
at energy E by the tj v (E) with the largest magnitude, 
which dictates convergence of the T matrix. This anal- 
ysis follows work in the early 1960's by Weinberg [48 1, 
so we call tj v a "Weinberg eigenvalue", although others 
have made similar treatments of the convergence of the 
Lippmann-Sch winger series. 

If we compare Eq. ( |48| ) to Eq. ( |47| ), we can express 
the convergence criterion for E < as saying that T(E) 
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Figure 36: Trajectories of Weinberg eigenvalues in the complex plane 
for several realistic NN potentials in two channels. The symbols go 
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diverges if there exists a bound state at E for V/rj v with 
\rjv\ ^ 1. Or, in other words, what is the largest r] v for 
which V/rjy supports a bound state at El We have con- 
vergence only for rj v < 1. This means we'll have two 
types of eigenvalues, because rj v could be negative ("re- 
pulsive") as well as positive ("attractive"). The negative 
eigenvalue corresponds to looking for bound states with 



a scaled "flipped" potential, as in Fig. 34 We see that 
the repulsive core becomes a deep attractive well, im- 
plying that we will have a large negative eigenvalue in 
these cases. But then we expect that RG evolution to a 
softened form, which eliminates the core, should result 
in decreased eigenvalues. 

These expectations hold in practice, as illustrated in 



Fig. 35 which shows the evolution of the largest posi- 
tive and negative Weinberg eigenvalues as a function of 
the V\ ow k cutoff A. (Very similar behavior is observed 
for the SRG with A replacing A.) In free space, the 
largest attractive eigenvalue is unity at all A's, corre- 
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sponding to the deuteron bound state. But the negative 
eigenvalue starts very large (because of the repulsive 
core) and drops dramatically as the potential is evolved, 
ending up well less than unity, indicating it is perturba- 
tive (but the positive eigenvalue still makes this channel 
nonperturbative). The situation is even more dramatic 
in the medium. The other curves in the figure show 
the result from considering the T matrix in the nuclear 
medium, where Pauli blocking effects are included in 
the intermediate states. We see a further reduction of 
the negative eigenvalue and now the positive eigenvalue 
is perturbative as well: the deuteron has dissolved! 

Starting from negative energies, we can follow the 
evolution of the largest eigenvalue into the complex 
plane as we increase E to positive values. Figures [36 



eigenvalues, with the symbols indicating energies of 0, 
25, 66, 100, and 150 MeV. The shaded region is the unit 
circle; the potential is perturbative for energies where 
both eigenvalues lie inside. The softening effect of the 



RG evolution is manifest in Fig. [37 and the eigenvalues 
provide a quantitative measure of the perturbativeness. 



It's also clear from Fig. 36 that at least the particular 
chiral potential used there is already quite soft. How- 



ever, Fig. [38] shows that significant additional soften- 
ing is possible with RG evolution. In all of these plots, 
the differences between the ^o channel and the 3 So- 
3 Do coupled channel stem from the latter having an ad- 
ditional source of nonperturbative behavior: the short- 
range tensor force. 

The increasing "perturbativeness" at finite density is 



and 37 show paths in the complex plane for Weinberg documented again in Fig. 39 At typical nuclear densi 



!— -,. 



0.5 



PQ 



-0.5 



R.J. Furnstahl / Nuclear Physics B Proceedings Supplement 00 (2013) l-\37\ 

150 



21 



S with Pauli blocking 



A = 4.0 fm 

A = 3.0fm 1 
A = 2.0fm ] 



J_ 



JL 



0.2 0.4 



0.6 0.8 
k F [fm" 1 ] 



1.2 1.4 



Figure 39: Dependence of the largest Weinberg eigenvalues on density 
for several evolved V\ ow k potentials |47 1. 



ties with 1 fm -1 < k$ < 1.3 fm -1 , both positive and neg- 
ative eigenvalues are small at the lowest A's. This im- 
plies that nuclear matter may actually be perturbative! 
(Note that at the Fermi surface, pairing as a nonpertur- 
bative phenomenon is revealed by \tj v \ > 1 |49|.) We 



can understand how this happens from Fig. 40 which 
shows the phase space available to two nucleons that 
scatter in the medium. Pauli blocking means they must 
go outside the two Fermi spheres, but the volume is in- 
creasingly restricted with decreasing A. In addition, the 
magnitudes of the matrix elements that scatter such par- 
ticles decrease as well 1 50 1 . 



A:IP/2±kl>£ F and Ikl <A 




F: IP/2 ± kl < k 



Figure 40: Overlapping Fermi spheres showing available phase space 
for two nucleons excited above the Fermi surface |50|. 



Perturbation theory in the particle-particle channel is 
shown in Fig. [41] for the high-resolution Argonne vig 
potential initially and after evolution by the Vi 0W £ RG 
to low resolution (A = 1.9 fm -1 ). Whereas many-body 
perturbation theory (MBPT) manifestly diverges for the 
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Figure 41: Many-body perturbation theory for symmetric nuclear 
matter up to third order in the particle-particle channel (33]. 



original potential, it converges for the low-momentum 
interaction at second order (at least in this channel; more 
complete many-body approximations must be studied to 
be more definite). It is also evident that there is no sat- 
uration. But adding a 3NF fit only to few-body prop- 



erties, as in Fig. [42] shows that the empirical saturation 
point can be reproduced with an uncertainty of about 2- 
3 MeV/particle. On-going work to improve this result 
includes the development of SRG evolution for the 3NF 
in momentum space l37l and of coupled cluster meth- 
ods for infinite matter with 3NF's to provide a high- 
order resummation of perturbation theory to test con- 
vergence. 

4.3. Universality from flow equations 

Another general aspect of RG flows known from the 
study of critical phenomena is the appearance of univer- 
sal behavior. In the application of RG to nuclear inter- 
actions, the universality we observe is that distinct ini- 
tial NN potentials that reproduce the experimental low- 
energy scattering phase shifts, are found to collapse to a 
single universal potential. We've already seen an indi- 
cation of this, but here we document it in more detail. 

We focus on the SRG, but very similar conclusions 
are found for Vi 0W £ evolution. In Figs. [43] and [44] S- 



wave N 3 LO chiral EFT potentials from Refs. |29] and 
lf5T1l are evolved with the SRG. Although the level of 
truncation is the same and the cutoffs approximately 
equal, the methods of regulating the potential differ 
(particularly for the two-pion exchange part). The result 
is very different looking initial interactions. This is not 
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a concern, because the NN potential is not an observ- 
able. We also observe there is significant off-diagonal 
strength coupling low and medium momenta in the ini- 
tial potentials. 

As the potentials are evolved, we see the characteris- 
tic driving toward the diagonal, with the diagonal width 
in k 2 given roughly by A 2 . At the end of the evolution 
shown, the interactions still look quite different at first 
glance. However, if we focus on the low-momentum 
region, where k 2 < 2 fm~ 2 , they appear much more sim- 
ilar. We can quantify this by taking a slice along the 
edge (i.e., V(k, 0) and along the diagonal (i.e., V(k, k)) 
and plot these quantities for these potentials and two ad- 
ditional ones. This is done in Fig. 45 We see a dra- 



matic collapse of the interaction between A = 5fm _1 
and A - 2fm _1 for the region of k below A (or maybe 
3/4 A). An open question under active investigation is 
whether evolved 3NF interactions will be universal. 

4.4. Operator evolution via flow equations 

We have focused almost exclusively on the evolution 
of the Hamiltonian, but an RG transformation will also 
change the operators associated with measurable quan- 
tities. It is essential to be able to start with operators 
consistent with the Hamiltonian and then to evolve them 
maintaining this consistency. The first step is a prime 
motivation for using EFT; we will assume that we have 
consistent initial operators in hand. The second step 
can be technically difficult, especially since we will in- 
evitably induce many -body operators as we evolve (for 
the same arguments as we made for the Hamiltonian). 
This is where the SRG is particularly advantageous, be- 
cause it is technically feasible to evolve operators along 
with the Hamiltonian. 



The SRG evolution with s (recall s = I /A 4 ) of any 
operator O is given by: 



O s = U s OU] , 



so O s evolves via 



dO s 

ds 



[[G S9 H s ],O s ] 



(50) 



(51) 



where we use the same G s to evolve the Hamiltonian 
and all other operators. While we can directly evolve 
any operator like this in parallel to the evolution of the 
Hamiltonian, in practice it is more efficient and numer- 
ically robust to either evolve the unitary transformation 
Ui itself: 



dU s 
ds 



= n s U s = [G s ,H s ]U s , 



(52) 



with initial value U s= o = 1 , or calculate it directly from 
the eigenvectors of H S = Q and H s : 



u s = Yj i<m*)><<m°)i • 



(53) 



Then any operator is directly evolved to the desired s 
by applying Eq. ( |50| ) as a matrix multiplication. The 
second method works well in practice. 

To simplify our study of operator evolution, we con- 
sider the simplest possible operator: the momentum 
number operator a q a q . In Fig. 46 we show the momen- 
tum distribution in the deuteron, i.e., < if/dWq^qWd > for 
two different realistic potentials, AVI 8 and CD-Bonn. 
As implied by the j-axis label, the momentum distri- 
bution is just the square of the deuteron wave function 
in momentum space. The results for the two potentials 
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Figure 43: SRG evolution of two chiral EFT potentials in the So channel |2 1. 
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Figure 44: SRG evolution of two chiral EFT potentials in the 3 Si channel 1 2 1. 
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Figure 45: SRG flow toward universality for several chiral EFT potentials |2|. 
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Figure 46: Momentum distribution of the deuteron for the AVI 8 po- 
tential and CD-Bonn potentials, compared to results from SRG poten- 
tials evolved to A = 2.0 fm -1 and 1.5 fm -1 . 
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Figure 47: Momentum distribution for an unevolved chiral N 3 LO 
potential compared to results from SRG potentials evolved to A = 
2.0 ml- 1 and 1.5 for 1 l52l. 



agree up to about 2fm _1 and then are different. If we 
evolve the AVI 8 potential and the momentum operator, 
then the matrix element < if/^O^i//^ > will be the same 
for any A and the curve will exactly reproduce the AVI 8 
deuteron momentum distribution. However, if we cal- 
culate < il/ A d \a\a q \il/ A d > for A = 2fm _1 and 1.5 fm -1 , that 
is, we use the evolved wave function but the unevolved 
operator, then we get the other curves. Besides directly 
illustrating that the high-g part of the momentum distri- 
bution is not an observable, since we can change it at 
will by unitary transformations, this manifests that the 
high-momentum part of the wave function is removed. 



The latter is potentially disturbing, because if the 
evolved operator is supposed to reproduce a high mo- 
mentum result when the evolved wave function has a 
vanishingly small component at that momentum, this 
may be because the operator is becoming pathological. 
To explore this further, we consider the momentum dis- 
at low (q = 0.34 fm -1 ) and moder- 



47 



tribution in Fig. 

ately high (q = 3.0 fm -1 ) values of the momentum [52]. 



(54) 



In Fig. 48 we plot the integrand of 

(i[/ d \(Uala q tf)\irt) , 

at each of these two q values. The full integral is the 
momentum distribution at those g's, so the plots tell 
us where the strength of the operator lies. For the 
low-momentum operator, there is little renormalization, 
but the nature of the high-momentum operator changes 
completely. Originally, the integral comes entirely from 
the region of q = 3.0 fm -1 , but the evolution of the oper- 
ator shifts its strength entirely to low momentum. This 
result is similar for other operators, such as electromag- 
netic form factors [52]. 

As we move to A > 3, the operator evolution and ex- 
traction process becomes more involved. A flowchart 

Imagine we ini- 



for the procedure is given in Fig. 49 



daily have a one-body operator and we want to evolve 
and then evaluate it in an A-particle nucleus. The dif- 
ficulty is that ft-body components are induced as we 
evolve and these must be separated out so we can cor- 
rectly embed them in the nucleus. In particular, to 
embed an ^z-body operator in an anti- symmetrized A- 
particle nucleus, we need an embedding factor of rM, so 
we need to isolate the components first. 

With the usual SRG generators, there is no evolu- 
tion of one-body operators. This is easiest to see in 
second quantization, where the commutators on the 
right side yield operators that have at least four cre- 
ation/destruction operators, meaning it is two- and 
higher-body. To isolate the two-body part, we first 
evolve the Hamiltonian in the two-particle basis to find 

(2) 

the unitary operator U s and use it to evolve our opera- 
tor. Then if we subtract the embedded one-body opera- 
tor, we will have our two-body part. The one- and two- 
body parts are then embedded in the 3-particle basis and 
subtracted from the evolved operator in that basis. And 
so on until we reach the desired level of truncation, at 
which point we embed in the A-particle basis and per- 
form the A-particle calculation. More details and exam- 
ples can be found about operator evolution in Ref. [52]. 

4.5. Computational aspects 

Before concluding this lecture, we'd like to make 
some brief comments on the computational aspects of 



R.J. Furnstahl / Nuclear Physics B Proceedings Supplement 00 (2013) l-\37\ 



25 



k'(fm" 1 ) k'(fm" 1 ) k'(fm" 1 ) k'(fm" 1 ) 

„0 123012301230123 



u 








* 


1 


_ 1 

£2 
" 3 


X =6.0 fm' 1 


X =3.0 fm' 1 


X =2.0 fm 


■ 1 


I 

X=1.5fm _1 



k'(fm" 1 ) k'(fm" 1 ) k'(fm" 1 ) k'(fm" 1 ) 

J) 123012301230123 



_ 1 

\ 2 
"3 



X =6.0 fm -1 



X =3.0 fm -1 



X =2.0 fm" 1 



X =1.5 fin" 1 



1 

0.5 



-0.5 

-1 

x10 3 

i 

0.5 


-0.5 
-1 



Figure 48: Integrand of momentum distribution operator |52|. 
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the calculations behind the figures we've seen. As noted 
earlier, the continuous momentum is discretized into a 
finite number of momentum points. The subsequent dis- 
cretization of integrals leads directly to matrices, and 
most of the manipulations are efficiently cast in this lan- 
guage. For example, momentum- space flow equations 
have integrals like: 



■/■ 



I(p, q)= \dk k 2 V(p, k)V(k, q) . 



(55) 



The usual choice of discretization is to use gaussian 
quadrature to accurately evaluate integrals with a mini- 
mum of points. (The size is not an issue for two-body 
operators, but becomes critical for higher-body opera- 
tors.) We introduce gaussian nodes and weights {k n , w n } 
with (n = 1, A/) to reduce integrals to finite sums: 

JdkmvJ^WnfiK). (56) 



(Note: these sets of nodes and weights are generally a 
combination of separate smaller rules over adjacent in- 
tervals with a total of N points.) Then I(p, q) — > I t j, 
where p = k t and q = kj, and 



// 



2> 



7 ? vv n Vin V n j 



Yj V inVn 



where 



Vij= ^ikiVijkj^j. 



(57) 



(58) 



This allows us to solve SRG equations and integral 
equations for phase shifts, Schrodinger equation in mo- 
mentum representation. In practice, N « 100 gauss 
points is adequate for accurate nucleon-nucleon partial 
waves. 

A computer code that carries out SRG evolution can 
be remarkable simple. Here is a possible pseudocode 
that is suitable: 



Such a code has been implemented in MATLAB, Math- 
ematica, Python, C++, and Fortran-90. While any basis 
can be used, so far only discretized momentum and har- 
monic oscillator bases have been implemented. Note 
that the same procedure (and even the same code in 
some cases) is relevant for a many-particle basis, but 
the number of differential equations will grow rapidly. 
For accurate results in the two-body evolution, 100 2 = 
10,000 matrix elements are needed, so there are that 
many differential equations. For an accurate 3NF evo- 
lution in a harmonic oscillator basis, at least 20 mil- 
lion coupled differential equations need to be solved. 
This sounds intimidating, but is well within reach of 
MATLAB (for example) on a machine with a moderate 
amount of memory. 

4.6. Summary points 

Chiral EFT establishes a hierarchy of many-body 
forces. Using flow equations to run to low resolution 
makes many-body calculations more perturbative and 
interactions flow to universal form (at least for NN; 
this is not yet established for the 3NF). Operators can 
be evolved consistently with interaction. Long-distance 
operators change very little, while short-distance oper- 
ators renormalize significantly, accompanied in some 
cases by a change in physics interpretation. The basic 
SRG flow equations in a partial wave momentum basis 
can be cast in a form that involves just matrix manipu- 
lations and the solution of ordinary first-order, coupled 
differential equations. 



5. Nuclear applications and open questions 

In this final lecture, we take a look at the in-medium 
similarity renormalization group (IM-SRG) and make a 
broad survey of nuclear applications. We conclude with 
a summary of open problems and new challenges. 



1. Set up basis (e.g., momentum grid with gaussian 
quadrature or HO wave functions with Af max ) 

2. Calculate (or input) the initial Hamiltonian and G s 
matrix elements (including any weight factors) 

3. Reshape the right side [[G S ,H S ],H S ] to a vector 
and pass it to a coupled differential equation solver 

4. Integrate V s to desired s (or A = s~ 1/4 ) 

5. Diagonalize H s with standard symmetric eigen- 
solver => energies and eigenvectors 

6. Form U = £i \ijr s )w* =0 \ from the eigenvectors 

7. Output or plot or calculate observables 



5.7. In-medium similarity renormalization group 

We start with a review of Hartree-Fock. The Hartree- 
Fock wave function is the best single Slater determinant 

|¥hf> = det{0;(x), i = 1 • • • A} , x = (r, a, r) (59) 

in the variational sense. The 0;(x) single-particle wave 
functions satisfy non-local Schrodinger equations: 



2M 7 



■/■ 



■i(x) + V R (x)4>i(x) + | dy V E (x, y)0Ky) = etfiix) 

(60) 



with direct 
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Wick's theorem can then be used to exactly write H as 
H = E + ^ f\2{a\a 2 } 



V H (x)= py£l0,(y)l 2 v(x,y) (61) 

J 7=1 



and exchange 



V E (i >y ) = -v(i,y)^/i)#)(y) 

7=1 



(62) 



potentials (54JI55). The direct and exchange potentials 



are shown diagrammatically in Fig. 50 with the right- 
most diagram an abbreviated form called a Hugenholtz 
diagram (56). We must solve self-consistently using oc- 
cupied orbitals for V H and Ve. Then Slater determinants 
from all orbitals form an A-body basis. 



o-o + (D 



Figure 50: Feynman and Hugenholtz diagrams for Hartree-Fock. 



12 



(2!) : 



l — J] (l2\T\34){a\ala 4 a 3 } 



1234 



(3!) : 



X — Yj (l23\T^\456){a\alala 6 a 5 a 4 },(64) 



123456 



where the zero-, one-, and two-body normal-ordered 
terms are given by 

E = <0|#|0> = 2>ii*i 

l 

+ j£<12|V|12>/!in 2 

Z 12 

+ ^Z <123|y(3)|123> ^^ 3 ' (65) 

' 123 

/i2 = r 12 + ^<i/|y|2/K- 

+ ^(10F (3) |2/jKn 7 -, (66) 



The in-medium SRG (IM-SRG) for nuclei, developed 
recently by Tsukiyama, Bogner, and Schwenk [57], ap- 
plies the flow equations in an A-body system using a 
different reference state than the vacuum. For exam- 
ple, we can choose the Hartree-Fock ground state as a 
reference state. The appealing consequence is that, un- 
like the free- space SRG evolution, the in-medium SRG 
can approximately evolve 3, ..., A-body operators using 
only two-body machinery. However, also in contrast to 
the free-space SRG, the in-medium evolution must be 
repeated for each nucleus or density. 

The key to the IM-SRG simplification is the use of 
normal-ordering with respect to a finite-density refer- 
ence state. That is, starting from the second-quantized 
Hamiltonian with two- and three-body interactions, 

H = Yj T na\a 2 + -r^ J] < 12 l^l 34 >^^3 



12 



(2!) : 



1234 



(3!) 2 



Y (l23\V (3) \456)a\alala 6 a 5 a 4 , (63) 



123456 



all operators are normal-ordered with respect to a finite- 
density Fermi vacuum |<D) (for example, the Hartree- 
Fock ground state or the non-interacting Fermi sea in 
nuclear matter), as opposed to the zero-particle vacuum. 



(12|r|34) = (12|V|34) + ^(12/|y (3) |34/)^, (67) 

i 

and rii - 0(sp-£i) denotes the sharp occupation numbers 
in the reference state, with Fermi level or Fermi energy 

By construction, the normal-ordered strings of cre- 
ation and annihilation operators obey (®\{a[ • • • <z„}|0) = 
0. It is evident from Eqs. ([65])— ([67]) that the coeffi- 
cients of the normal-ordered zero-, one-, and two-body 
terms include contributions from the three-body interac- 
tion V (3) through sums over the occupied single-particle 
states in the reference state |<D). Therefore, truncat- 
ing the in-medium SRG equations to two-body normal- 
ordered operators will (approximately) evolve induced 
three- and higher-body interactions through the density- 
dependent coefficients of the zero-, one-, and two-body 



operators in Eq. ( 64 ). 

The in-medium SRG flow equations at the normal- 
ordered two-body level are obtained by evaluating 
dH/ds = [77, H] with the normal-ordered Hamiltonian 
H = £" +/+r and the SRG generator n = n lb +n 2b (with 
one- and two-body terms) and neglecting three- and 
higher-body normal-ordered terms. For infinite mat- 
ter, a natural generator choice is n = [/, T] in analogy 
with the free-space SRG. In this case, the explicit form 
of the SRG equations simplifies because n lb = and 
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fij = ft 8ij. This leads to 



~di ~ 2^ (fl2 ~ M 



)|r 123 4l nin 2 n 3 n 4 , (68) 



1234 



—j- = / j CAl - /23) 1^4123 1 



234 



x (n^n^n^ + n^n^n^ , 



(69) 



— ~, - -(J12 - /34J 1 1234 

as 



+ ~Z 2_j (/l2 + /34 _ 2fabWl2abT a b34 

ab 

x(l-n a -n b ) 
+ £j(n a - n b ) 

ab 
X \^alb3^b2a4[ifal ~ fbl) ~ (fb2 ~ fa*)] 
~ r a2 b3r b la4[(fa2 ~ fbs) ~ (fbl ~ fa*)]} , (70) 

where the single-particle indices refer to momentum 
states and include spin and isospin labels. 

While the in-medium SRG equations are of second 
order in the interactions, the flow equations build up 
non-perturbative physics via the successive interference 
between the particle-particle and the two particle-hole 
channels in the SRG equation for T, Eq. ( [70} , and 
between the two-particle-one-hole and two-hole-one- 
particle channels for /, Eq. ( [69] ). In terms of diagrams, 
one can imagine iterating the SRG equations in incre- 
ments of Ss. At each additional increment Ss, the in- 
teractions from the previous step are inserted back into 
the right side of the SRG equations. Iterating this pro- 
cedure, one sees that the SRG accumulates complicated 
particle-particle and particle-hole correlations to all or- 
ders. 

With the choice of generator r] = [f,T], the Hamil- 
tonian is driven towards the diagonal. This means 
that Hartree-Fock becomes increasingly dominant with 
the off-diagonal T matrix elements being driven to 
zero. As with the free-space SRG, it is convenient for 
momentum- space evolution to switch to the flow vari- 
able A = s~ 1/4 , which is a measure of the resulting band- 
diagonal width of T. In the limit A — > 0, Hartree-Fock 
becomes exact for the evolved Hamiltonian; the zero- 
body term, Eq, approaches the interacting ground- state 
energy, / approaches fully dressed single-particle ener- 
gies, and the remaining diagonal matrix elements of T 
approach a generalization of the quasiparticle interac- 
tion in Landau's theory of Fermi liquids 1 1 ]. 

An approximate solution of the Eq flow equation for 
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as evolved in the IM-SRG ED. 
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1.35 fm as 



symmetric nuclear matter and neutron matter as a func- 



tion of A is shown in Figs. 51 and 52 for two differ- 
ent Fermi momenta k F (corresponding to different den- 
sities). As expected, the in-medium SRG drives the 
Hamiltonian to a form where Hartree-Fock becomes 
exact in the limit A — » 0. In contrast to the lad- 
der approximation based on NN-only SRG interactions 
evolved in free space, the same many-body calcula- 
tion using interactions evolved with the in-medium SRG 
at the two-body level gives energies that are approxi- 
mately independent of A. This indicates that truncations 
based on normal-ordering efficiently include the dom- 
inant induced many-body interactions via the density- 
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Figure 53: "Off-diagonal" terms (e.g., 2p2h sectors) driven to zero as s increases, decoupling them from the Hartree-Fock reference state, which 
becomes exact as s — » oo |57l . 
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Figure 54: IM-SRG flow for the energy of the alpha particle (57l . 
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dependent zero-, one-, and two-body normal-ordered 
terms. 

In a similar manner, the in-medium SRG can be used 



as an ab initio method for finite nuclei. Figure 53 shows 
the off-diagonal normal-ordered two-body matrix ele- 
ments r being driven to zero with the IM-SRG evolu- 



tion. Figure 54 shows the evolution of the ground-state 
energy of 4 He [57]. As the flow parameter s increases, 
the Eq flow and second-order (in T) many-body per- 
turbation theory contributions approach each other, as 
was the case for the infinite matter results in Fig. |5T 
In addition, the convergence behavior with increasing 
harmonic-oscillator spaces in Fig. 55 for 4 He and 40 Ca 
is very promising. Based on these calculations, the in- 
medium SRG truncated at the normal-ordered two-body 
level appears to give accuracies comparable to coupled- 
cluster calculations truncated at the singles and doubles 
(CCSD) level. Finally, we note that the in-medium SRG 
is a promising method for non-perturbative calculations 
of valence shell-model effective interactions and opera- 
tors. 

5.2. Implications ofRGfor nuclear calculations 

There are many on-going and potential applications 
of low resolution methods for calculations of nuclear 
structure and reactions. In many cases RG techniques 
are used explicitly but there are also examples of low 
resolution being achieved by other means. Here we'll 
survey some of both types. This will be far from a com- 
prehensive list because it focuses largely on develop- 
ments associated entirely or in part with a project called 
UNEDF, which stands for Universal Nuclear Energy 
Density Functional. 

UNEDF is a collaboration of more than fifty physi- 
cists, applied mathematicians, and computer scientists 
in the United States plus many international collabora- 
tors, funded through the U.S. Department of Energy's 
SciDAC program. The long-term vision of the project 
is to arrive at a comprehensive and quantitative descrip- 
tion of nuclei and their reactions. The focused mis- 
sion is to construct, optimize, validate, and apply en- 
ergy density functionals for structure and reactions, but 
to carry out this mission the team has developed many 
crosscutting physics collaborations where none existed 
previously between the main physics areas: ab initio 
structure, ab initio functionals, DFT applications, DFT 
extensions, reactions. These interconnections are indi- 
cated schematically in the UNEDF strategy diagram in 
Fig. [56] This type of large-scale collaboration repre- 



sents a transformation in how low-energy nuclear the- 
ory is done. UNEDF has been very productive, with 



over 200 publications to date, including 1 1 Physical Re- 
view Letters and a Science article in the 2011 calendar 
year alone. Further background, references, and scien- 
tific highlights can be found at unedf.org, the project 
website. 

One of the important tools for nuclear structure en- 
abled by low-resolution interactions is the diagonaliza- 
tion of enormous but very sparse Hamiltonian matrices, 
usually in a harmonic oscillator basis to permit center- 
of-mass effects to be excluded. This is referred to as no- 
core full configuration (NCFC) or no-core shell model 
(NCSM), depending on the context. Two recent ex- 
amples of what is enabled are represented in Fig. 57 
On the left are Gamov-Teller matrix elements from a 
large-scale calculation of Carbon- 14 using a soft chiral 
EFT potential [58]. They highlight the critical role of 
the 3NF in suppressing the beta decay rate, which ex- 
plains the anomalously long lifetime of 14 C (which is 
used to great advantage for dating artifacts!). On the 
right is the low-lying spectrum of Fluorine- 14, which 
is unstable to proton decay. This spectrum was pre- 
dicted in advance of the experimental measurements 
(not a common occurrence until now!), which meant 
solving a Hamiltonian matrix of dimension 2 billion us- 
ing 30,000 cores with a soft interaction (derived from 
inverse scattering rather than RG, but with similar char- 
acteristics) [59, 60]. The predictions and measurement 
agree within the combined experimental and theoretical 
uncertainties. These calculations would not be possible 
with the potentials of Fig. [3] RG- softened interactions 
will allow many more of these confrontations of exper- 
iment with theory in the future. 

One of the principal aims of the UNEDF project is 
to calculate reliable reaction cross sections for astro- 
physics, nuclear energy, and national security, for which 
extensions of standard phenomenology is insufficient. 
The interplay of structure and reactions is essential for 
a successful description of exotic nuclei as well. Such 
interplay is characteristic of the ab initio no-core shell 
model/resonating-group method (NCSM/RGM), which 
treats bound and scattering states within a unified frame- 
work using fundamental interactions between all nucle- 
ons. A quantitative proof-of-principle calculation of this 
approach is shown in Fig. [58] [61 , 62]. A wide range 
of applications is now possible including 3 H(d,ft) 4 He 
fusion and the 7 Be(/?, y) 8 B reaction important for solar 
neutrino physics, and many more to come. In Fig. [59 
the first-ever ab-initio calculation of the 7 Be(/?, y) 8 B as- 
trophysical S-factor is shown l63ll . This calculation uses 
NCSM/RGM with an N 3 LO NN interaction evolved by 
the SRG to a fine-tuned value of A = 1.86 fm" 1 . The 
ab initio theory predicts both the normalization and 
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Figure 58: NCSM/RGM calculation of neutron scattering from the 
alpha particle using SRG-evolved interactions, compared with exper- 
imental measurements 1 61 62 1 . 
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the shape of 5 17. This very promising techniques has 
many additional applications. The inclusion of an SRG- 
evolved 3NF is planned for the near future. 

UNEDF members have demonstrated that the pow- 
erful coupled-cluster (CC) ab initio method, which is a 
workhorse in quantum chemistry, can be used to accu- 
rately calculate closed-shell medium-mass nuclei such 
as 40 C, 48 C, 56 Ni with chiral EFT two-body interactions 



or the RG-softened versions [64, 65], as well as proton 
halo nuclei like 19 F [66]. A proof-of-principle conver- 
gence curve (ground-state energy versus the size of the 
orbital space) for 56 Ni with an SRG-evolved potential is 



shown in Fig. 60 The CC formalism has been extended 
to include NNN forces and their inclusion in calcula- 
tions of the heavier nuclei will break new barriers. 
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Figure 60: Convergence of CCSD for Ni-56 evolved N3LO to A = 
2.5 for 1 . 



The in-medium SRG diagonalization of clo sed- s hell 
nuclei such as 40 C (57), discussed in Section 5.1 



is a 

complementary approach to CC and is one of several 
advances in our understanding of the phenomenological 
nuclear shell model enabled by softened potentials. An- 
other is the direct use of MBPT to examine the effect 
of the 3NF on the location of the neutron dripline: the 
limits of nuclear existence where an added neutron is no 
longer bound — it "drips" away. The new physics is indi- 
cated schematically in Fig. [61] It's been established ex- 



perimentally that as you add neutrons to stable oxygen- 
16, the neutrons stay bound until 24 0. But adding one 
more proton to get fluorine extends the dripline all the 
way to 31 F. This result is not predicted by previous mi- 
croscopic calculations using NN interactions, because 
the single-particle neutron energy levels that get filled 
are predicted to be bound (leftmost panel in Fig. |62| ), 
leading to 28 as the calculated dripline. The phe- 
nomenological shell model, in which matrix elements of 
the Hamiltonian are fit to nearby nuclei, has a very dif- 
ferent pattern (e.g., compare the d3/2 single-particle en- 



16 



O core 




Figure 61 : Influence of a three-body force on valence neutrons in oxy- 
gen isotopes. 
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ergies in the middle and left panels of Fig. 62 ) and pre- 
dicts the correct dripline. However, recent calculations 
with a Viowfc RG force and fitted 3NF yield the right 
panel of Fig. [62] When the 3NF effect is added 



the interaction of valence neutrons with a core neutron, 
as in Fig. [61] is repulsive, pushing up the d3/2 level so 
that the dripline is at 24 0. 

The apparent success of MBPT with low-momentum 
potentials has been tested by Roth and collabora- 
tors [68], who have done calculations of closed shell 
nuclei across the mass table in second-order perturba- 
tion theory (first-order is Hartree-Fock). Results for 
SRG evolved interactions from an initial NN chiral EFT 
and including the induced 3NF show excellent indepen- 
dence of the flow parameter A (see Fig. [63]) and both the 
energies and radii are in good agreement with coupled 
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SRG-evolved NN interactions, including the induced 3NF |68|. Sev- 
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cluster results (64JI63. However, adding a 3NF leads to 
large X dependence. Very recent results from Roth et al. 
that use importance truncated NCSM as well as coupled 
cluster calculations show that the long-range 3NF is the 
source to apparent large 4NF contributions for oxygen 
and heavier nuclei, causing a strong dependence on the 
flow parameter. However, by using a lower cutoff for 
the initial 3NF, remarkable agreement with experimen- 
tal binding energies is achieved with fits only to few- 
body properties (19] [41]]. Work is in progress to identify 
SRG generators to better control the RG evolution of the 
initial 3NF. 

We have already seen the convergence of MBPT for 
symmetric infinite nuclear matter. Perturbation theory 
is even more controlled for pure neutron matter, as il- 
lustrated in Figs. [64] and [65] where Vi ow & RG-evolved 



interactions and fitted 3NF's are used in calculations of 
the neutron matter energy per nucleon as a function of 
the density 1 69 ]. The cutoff dependence of the result is 
used to estimate the many-body uncertainty. Figure 64 
shows that the 3NF contribution is important (compare 
to the NN-only curves) but that the dominant theoret- 
ical uncertainty is the value of the coupling constants 
for the long-range part of the N 2 LO chiral EFT 3NF. In 



Fig. 65 comparisons with non-perturbative calculations 
demonstrate the consistency of the much easier MBPT 
calculations. The results for the neutron matter equation 
of state have been used by Hebeler et al. to provide tight 
constraints on neutron star masses and radii ITOll . 

The MBPT results for infinite matter are also valu- 
able input for work to develop a microscopic nuclear 
energy density functional (EDF). A key tool to incorpo- 
rate microscopic input into nuclear EDF's is the den- 
sity matrix expansion (DME) originally proposed by 
Negele and Vautherin, which has been revived and im- 
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compared to other calculations 1 69 |. 



proved in the UNEDF project. The DME provides a 
route to an EDF based on microscopic nuclear interac- 
tions through a quasi-local expansion of the energy in 
terms of various local densities and currents, including 
resummations that can treat long-range one- and two- 
pion exchange interactions given by chiral EFT. With 
sufficiently soft microscopic interactions, many-body 
perturbation theory (MBPT) for nuclei is a quantitative 
framework for implementing the DME. The formal de- 
velopment of DME with MBPT is on-going, but there 
are already hybrid formulations between purely ab ini- 




Figure 66: Schematic representation of interacting neutrons in a (the- 
oretical) trapping potential. 



tio and phenomenological functionals [71], which allow 
improvements to be made while more systematic func- 
tionals are developed. 

Several DME implementations strategies have been 
developed, with the first tests recently made against 
ab-initio calculations using a semi-realistic interaction 
(Minnesota) in trapped neutron drops (73. Neutron 
drops are a powerful theoretical laboratory for improv- 
ing existing nuclear energy functionals, with particu- 
lar value in providing microscopic input needed for 
neutron-rich nuclei, where there are fewer constraints 
from experiment. The necessity of an external potential 
(because the untrapped system is unbound, with pos- 
itive pressure) is turned into a virtue by allowing ex- 



ternal control over the environment (see Fig. 66 for a 
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schematic of the trapped neutrons). Density functional 
theory, which provides the theoretical underpinning for 
the microscopic EDF's, dictates that the same functional 
applies for any external potential, which can therefore 
be varied to probe and isolate different aspects of the 



EDF. Results summarized in Fig. 67 show promising 
agreement with the DME functionals and the (essen- 
tially exact) ab initio results using NCFC. Many more 
developments in this line will be forthcoming. 

5.3. Summary and survey of open problems 

In these lectures, we've made a whirlwind tool of 
atomic nuclei at low resolution. With the renormaliza- 
tion group (RG), the strategy has been to lower the reso- 
lution and track dependence on it. We've seen how high 
resolution leads to coupling of low momenta to high 
momenta, which hinders solutions of the many-body 
problem for low-energy properties. With RG evolution, 
correlations in wave functions are reduced dramatically, 
leading to faster convergence of many-body methods. A 
consequence is that non-local potentials and many-body 
operators are induced, so these must be accommodated. 

Flow equations (SRG) achieve low resolution by de- 
coupling. This can be in the form of band or block diag- 
onalization of the Hamiltonian matrix. The flow equa- 
tions implement a series of unitary transformations, in 
which observables (measurable quantities) are not al- 
tered but the physics interpretation can change! In the 
nuclear case, the usual plan is to evolve until few-body 
forces start to grow rapidly, or to use an in-medium ver- 
sion of the SRG. 

With the RG, cutoff dependence becomes a new tool 
in low-energy nuclear physics. The basic idea is that, 
in principle, observables should be unchanged with RG 
evolution. In practice, there are approximations in the 
RG implementation and in calculating nuclear observ- 
ables. These come from truncation or approximation of 
"induced" many-body forces/operators and from many- 
body approximations. For nuclei there can be dramatic 
changes even with apparently small changes in the reso- 
lution scale. We can use these changes as diagnostics of 
approximations and to estimate theoretical errors. Some 
specific applications of cutoff dependence include: 

• using cutoff dependence at different orders in an 
EFT expansion, which carries over to the corre- 
sponding RG-evolved interactions; 

• using the running of ground-state energies with 
cutoff in few-body systems to estimate errors and 
identify correlations (e.g., the Tjon line); 



• in nuclear matter calculations, validating MBPT 
convergence and setting lower bounds on the er- 
rors from uncertainties in many-body interactions; 

• in calculations of finite nuclei, diagnosing missing 
many -body forces; 

• identifying and characterizing scheme-dependent 
observables, such as spectroscopic factors. 

The possibilities have really only begun to be explored. 

We have seen glimpses of the many promising ap- 
plications of RG methods to nuclei. Configuration in- 
teraction and coupled cluster approaches using softened 
interactions converge faster, opening up new possibil- 
ities and allowing the limits of computational feasibil- 
ity to be extended. Ground-breaking ab-initio reaction 
calculations are now possible. Applications of low- 
momentum interactions to microscopic shell model cal- 
culations bring new understanding to phenomenological 
results, highlighting the role of three-body forces. Be- 
cause many-body perturbation theory (MBPT) is feasi- 
ble with the evolved interactions, the door is opened to 
constructive nuclear density functional theory. 

There are also many open questions and difficult 
problems in applying RG to low-energy nuclear physics. 
Here is a subset: 

• Power counting for evolved many-body operators. 
That is, how do we anticipate the size of contri- 
butions from induced many-body interactions and 
other operators? This is essential if we are to 
have reliable estimates of theoretical errors, be- 
cause truncations are unavoidable. We need both 
analytic estimates to guide us as well as more ex- 
tensive numerical tests. Many of the same issues 
apply to chiral EFT; can the additional information 
available from SRG flow parameter dependence 
help with analyzing or even constructing EFT's? 

• Only a few possibilities for SRG generators have 
been considered for nuclear systems. Can other 
choices for the SRG G s operator help to control the 
growth of many-body forces? Can convergence be 
improved in the harmonic oscillator basis, which 
is limited by an infrared cutoff as well as an ultra- 
violet cutoff? Can a generator be found to drive 
non-local potentials to local form, so they can be 
used with quantum Monte Carlo methods? Or can 
the SRG equations be formulated to directly pro- 
duce a local projection and a perturbative residual 
interaction? 

• An apparent close connection between the block- 
diagonal generator SRG and the "standard" Vi ow fc 
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RG has been established empirically, but a formal 
demonstration of the connection and its limits has 
not been made. 

• What other bases for SRG evolution would be ad- 
vantageous? The need for a momentum- space im- 
plementation for evolution in the A = 3 space and 
beyond is foremost. Beyond providing necessary 
checks of evolution in the harmonic oscillator ba- 
sis, the evolved interactions in this form could be 
directly applied to test MBPT in infinite matter and 
to test nuclear scaling. Another possibility is to use 
hyperspherical coordinates, which combine the ad- 
vantage of a discrete basis with better asymptotic 
behavior (and which would be useful for visualiza- 
tion of many-body forces). 

• There are many open questions and problems in- 
volving operators. These include formal issues 
such as the scaling of many-body operators and 
technical issues such as how to handle boosts of 
operators that are not galilean invariant. And there 
are simply many applications that are yet to be 
made (e.g., electroweak processes). 

• The flow to universal form exhibited by two-body 
interactions has been clear from the beginning of 
RG applications to nuclei, but whether this same 
behavior is expected for many -body interactions or 
for other operators is still open. 

• How can we use more of the power of the RG? 
There is no shortage of opportunities and challenges ! 
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